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In the standard approach to the neutrino oscillations a unitary relation 
among weak and mass eigenstates of light neutrinos is imposed. However, 
O i| in many extensions of the SM left-handed, active neutrinos mix with ad- 

1} . ditional heavy neutrino states. Consequences of this additional mixing, 

driven by experimental constraints, on the neutrino oscillations are consid- 
ered. 

•T— I 

X. 

1. Introduction 

At present, 3 light neutrinos, with masses at the eV or sub-eV scale 
[1] are known to exist. However, much heavier neutrino states (m^r > 
0(Mz/2)) are not excluded [2]. These, due to kinematical reasons do not 
contribute directly to the weak neutrino states which can undergo neutrino 
oscillations. They influence, however, neutrino oscillations since they modify 
the neutrino mixing matrix U. Let TJ V be the full neutrino mixing matrix, 
then the matrix U of dimension 3x3 constitutes the mixing submatrix of 
light neutrino states (f e , u^, u T <-> u\, v^, z/3 transitions) 
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u »=(v' u>)- W 

The submatrix V (of dimension (tir x 3)) is responsible for the mixing of 
light neutrinos with hr heavy states ( v e , v^, u T <-> 1/4, . . . , v nR -3 transitions). 
The submatrix U' (of dimension {ur x tir)) is responsible for mixing among 
heavy states. In the conventional see-saw mechanism >> m u , where 
itin and m v are masses of heavy and light states respectively, the elements 
of V are very small and U becomes unitary. This simply means that heavy 
neutrino states do not modify mixings among light neutrino states. From 
the theoretical point of view, V must not necesserily be negligible [3]. We 
will use the experimental data to constrain V [4], and more precisely the 
combination ( VVM 1 (a,(3 = {e,fi, r}). 

V /a/3 

From the unitarity of U v we infer that 

The aim of this paper is to examine the effect of this modification of 
unitarity of V on neutrino oscillations. The subject is not new [6] 2 . Never- 
theless, some issues, especially connected with CP-violation effects have not 
yet been discussed. CP violation effects in the unitary neutrino oscillations 
case are known to be very fragile. If any element of the unitary U matrix 
(e.g. Ue^) is small then the effect of CP violation will be small either. And 
in fact, U e 3 (see Eq. 6) is known to be very small if not zero. Besides, the 
CP phase sin 5 must be substantial. Finally, the CP violating effects vanish 
with decreasing 5rriQ. For 5 = ^, <5toq given by LMA MSW solution and 
U e 3 > 0, the CP effects can be detectable [9], but even then it may hap- 
pen that matter effects will mimic (or screen) the CP violation [10]. We 
show that the nonunitarity of U can be responsible for similar effects. If 
CP-violation effects were detected with a strength larger than predicted by 
the unitary neutrino mixing approach, then heavy neutrino mixing could be 
held responsible for this effect. In the contrary case, some better bounds on 
the ( VV"n factors could be found. 

V /a/3 

In this paper we focus on neutrino oscillations in vacuum. 



1 The elements of the V matrix can also be investigated, e.g. in heavy neutrino pro- 
duction processes [5]. 

2 Lately, effects of a non-unitary mixing matrix U have been considered in [7] in a 
different context where new leptonic interactions have been included. 
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2. Neutrino oscillations in the presence of heavy neutrino states 



In the standard neutrino oscillation theory of three flavours we start 
with neutrino weak eigenstates v a = (i/ e , v^v T ) as a combination of three 
mass eigenstates i/j = (^1,^2,^3) 



3 

£ 

i=i 



(3) 



The form of the matrix U can be obtained using subsequent rotations 
around the axes spanned by massive neutrino states mi, m2, 



U — i?23-Rl3-Rl2- 



(4) 



.Ry's represent rotations in the i-j plane by 0jj angle with additional 
phases, e.g. (c i2 = cos9 12 , S12 = sin Q 12 e tSl2 ) 



R 



12 




(5) 



Taking 812 = (J23 = (two of the three complex phases do not influence 
the oscillation probability [11]) we obtain the classical parametrization of 
the U matrix [12] (S 13 = 5) 



U 




( 

V 



C12C13 

-C23S12 - Sl3«23Cl2e 
Sl2«23 - Sl 3 C 2 3Cl2e 



16 



C13S12 
C12C23 - Sl2-S23«13e i<5 
-S23C12 - Sl2C23Sl3e* 



si3e 

C13S23 I (6) 
C23C13 



Let us now include effects of the matrix V to the matrix U (Eqs. 1,2). 
We will do it by introducing three new parameters e^, i=l,2,3 which are 
connected directly to the elements of the matrix V in the case of the 4x4 
matrix Eq. 1. 

The general 4x4 matrix Eq. 1 can be parametrized by 6 rotation angles 
(and 6 phases) in the following way 

U v = R34R24R14R23R13R12, (7) 
where the rotations take place in the 4 dimensional space spanned by 
four massive neutrino states, e.g. [13], [14] (si 2 = sin Q\2e %Sl2 ) 



R12 



( 


C12 


6 12 





°\ 




-S12 


C12 
















1 





V 











1/ 



(8) 
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Let us note that Eq. 7 differs from Eq. 4 by three additional rotations 
-R34-R24-R14 in the plane to which the additional fourth neutrino state be- 
longs. When the fourth state is much heavier than the light states, the 
rotation angles are small. Let us take | S14 1 = |ei| << 1, | S24 1 = |e 2 | « 1 
and I S34 1 = I £3 1 << 1. Then we can expand Eq. 7 to get 



/ 



U v = 



Ufa 



\ 



V 9 fa) 



1 



ei 

C3 



ki| 2 + |6 2 | 2 + |e 3 | 2 ) / 



(9) 



In the limit e» -> Ufa) -> U (Eq. 6) and gfa) -> 0. Ufa) is the 
desired matrix which we will use in the neutrino oscillation formula instead 
of the U matrix in Eq. 3. We will not show the manifest form of Ufa) as it 
is straithorward but space consuming. Our parametrization through (com- 
plex) e factors holds in the general case of n heavy states and can be easily 
connected to the quantities which are usualy constrained by experimental 
data, e.g. 



x; \ v « 

-heavy 



hr < 0.0054, 



(10) 



and similarly, 



VV*) I = |eie 2 | < 10" 



(vv-t) 



/IT 



|e 2 e 3 | < 10" 



(11) 
(12) 



The very strict constraint, Eq. 11 comes from the lack of the fi — > e^f 
decay [4, 6, 15]. Constraints, Eqs. 10,12 are consequences of global fits to 
experimental data [4, 6] (e.g. lepton universality, invisible Z decay, CKM 

unitarity). There is also a constraint on \Vei\ 2 /Mi coming from the 

i=heavy 

neutrinoless double beta decay [8]. In our approach we do not have to use 
any information on the heavy neutrino mass spectrum. We just assume that 
the masses are above 100 GeV. The constraint from (f3(3)o u is then fulfilled. 
With the parametrization Eq. 9 it is straightforward to write the modified 
neutrino oscillation probability 



NlN} 



\ ) at 
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- 4 Rfp sin2 A ab ~ 8/^ sin A 2 i sin A 31 sin A 32 

a>b 



4 1 jsin2A3i+^sin2A 3 2 



i (2) 



(13) 



where 



A ab = 1.27Sm 2 ab [eV] 



L[km] 
E[GeV] 



5m 2 ab = m 2 a -m 2 b . 



(14) 



R^q are modified definitions taken from the standard, unitary approach 



R 



ab 

a(3 



fab 
l a/3 



Re 
Im 



W. 



a 3 



(15) 

W$,(ei) = U aa U ph U* ab U* pa . (17) 



W. 



N a (/3) are factors which normalize the three light neutrino states to 1 



N 2 = - 

1 - (Wt) 



The last row in Eq. 13 with 



Im 



all 



(18) 



(19) 



deserves an extra comment. Its appearence is a consequence of the mod- 
ification of the Jarlskog factors, which for unitary U fullfil the following 
relations 



jab 

I a/3 



jba 



rab jba 

_i /3a — 1 /3a- 



(20) 



When U is not unitary, Eq. 2 leads to 



712 



721 
l a/3 



723 
L af3 



732 
l a/3 



731 
l af3 



Im [U* a2 Up 2 (vV^) afj 

im (u* al u pi (V0 



(21) 
(22) 
(23) 
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and therefore to the last term in Eq. 13. 

As discussed in [6], the effects of the normalization factors will be diffi- 
cult to observe in experiments. Here we will focus on the influence of the 
additional neutrino mixing of light neutrinos represented by the e^'s on the 
CP violating effects. The novelty here is the appearance of the third line in 
Eq. 13. This term is not very sensitive to A21 when it is small. Therefore 
CP violation can occur even if <5mf 2 = 0. However, CP violation is now pos- 
sible with two neutrino oscillations. In addition, the CP effect with three 
neutrino flavours, contrary to unitary oscillations, can be substantial even 
if one of the elements of the mixing matrix is very small. 

3. CP-violating effects in neutrino oscillations 

CP violating effects can be seen in appearence experiments. Let us 
consider the following standard quantities 



A C p{a;(J) = — — , 

P(v a -> up) + P(u a -> up) 

P(u a -> up) - P{yp -> u a ) 



A T (a;P) 



P(u a -> up) + P{yp -»■ u a ) 



In vacuum Acp = At- The same is true in the case of the new Eq. 13 
when a nonunitary matrix U is present. 

We assume the following values of the standard parameters 



8m 2 21 = 5 • 1CT 5 eV 2 , 
5m 2 31 = 3 • 1CT 3 eV 2 , 
e i2 ^ 35°, 6 23 ^ 40°, 0i3 ~ 5°, (24) 
5 ~ ±90°. 

These values are consistent with CHOOZ [16], the LMA MSW solution 
of the solar neutrino problem [17] and the Superkamiokande data [18]. For 
the nonstandard parameters we take 



ei| = 0.001 

e 2 | = 0.1 (25) 

€ 3 1 = 0.1 



which are consistent with Eqs. 10,11,12. 
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Figs. 1,2 show the results for two cases of Acp(e; fi) and Acp(fi; r), and 
long-baseline (L=732 km, e.g MINOS) or short-baseline (L=250 km, e.g. 
K2K) experiments. The neutrino energies are chosen to be between 2 GeV 
and 30 GeV. We can see that the effects of the nonstandard heavy neutrino 
mixings can be quite large, even much bigger than in the unitary approach 
when 6i = 0. 

In Figs. 3,4 the results are given for genuine CP effects of NS sector 
when some of the e^'s are chosen to be complex and 5 = 0. 

Three lessons can be learned from these numerical results. First of all 
short baseline experiments are sensitive to the NS sector. Some improve- 
ments of the constraints Eqs. 10,11,12 are possible in this case when no 
signal for Acp is found. Second, the NS effects connected to the complexity 
of €i can mimic SM effects of 5. Third, cancellations between the SM and 
NS effects can appear. 

We would like to finish with a somehow academic example of what more 
do 'nonorthogonaP neutrino states mean. 

When a unitary U is used in the description of neutrino oscillations, the 
following relation holds 

E P ^ = 1 ' ( 26 ) 

a 

e.g.: 

Pee ~)~ Pefi "f" Per = 1- 

It simply means that the number of emitted neutrinos of the given 
flavour (e.g. e) will be the same as the number of final neutrinos of any 
type. However, for a nonunitary U this relation is not fullfil. Let us see it 
in a simple case of two flavours, when U is defined as (©2 = ©i + e) 

U=( COS@ ( \ Shl f}). (27) 
\ — sin ©2 cos ©2 / 

In this case we get 

Pec = Pee + -P eM = 1 + 4e sin 2 A 2 i sin ©1 cos ©1 cos 2©i + C(e 2 ](28) 

ct=e,fj, 

P ^ = p v* + p w = l-4esin 2 A 2 i sin ©! cos ©i cos 2©! + 0(e 2 ()29) 

a=e,/j, 

We can see that the sum can be either larger or smaller than 1. A similar 
result holds for a 3 dimensional U. 
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Fig. 1. The Acp(n;e) asymmetry as function of neutrino energy. The label 'NS' 
means that Eqs. 24,25 are taken into account. The label 'SM' means that the 
values of neutrino parameters as given in Eq. 24 and a = have been taken. 
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Fig. 2. The Acp(h;t) asymmetry as function of neutrino energy. The label 'NS' 
means that Eqs. 24,25 are taken into account. The label 'SM' means that the 
values of neutrino parameters as given in Eq. 24 and a = have been taken. 
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L=250 km, pure complex e, or£. 
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Fig. 3. The Acp(fi; e) asymmetry generated by the NS sector. The results are for 
the parameters Eq. 24 but with S = 0. ei = 0.001 • (1 or i), e 2 = 0.1 • (i or 1), 
£3 = 0.1. This choice is consistent with Eqs. 10-12. 
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Fig. 4. The Acp(p;r) asymmetry generated by the NS sector. The results are for 
the parameters Eq. 24 but with 6 = 0. e 2 ( r 3) = 0.1 • i, ei = 10~ 4 . This choice is 
consistent with Eqs. 10-12. 



